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Abstract. We consider different variants of factorization of a 2x2 matrix Schroedin- 
ger/Pauli operator in two spatial dimensions. They allow to relate its spectrum to the 
sum of spectra of two scalar Schrodinger operators, in a manner similar to one-dimensional 
'sj- ' Darboux transformations. We consider both the case when such factorization is reduced 

' to the ordinary 2-dimensional SUSY QM quasifactorization and a more general case which 

. involves covariant derivatives. The admissible classes of electromagnetic fields are described 

I and some illustrative examples are given. 

o 

O^; 1. Introduction 

In the ordinary one- dimensional supersymmetric quantum mechanics (SUSY QM) 
^ ! [1], a pair of Hamiltonians (superpartners) is factorized as: 

= -d^ + {dx?-d\ = s+s- d = 4- 

ax 

i/W = -d^ + {^dxf + d\ = S-S+ (1) 

where 

S^^^d + dx, 5+ = (5-)t. (2) 

The arbitrary function x{^) is called a superpotential. The (not necessarily normal- 
izable) solution of the Schrodinger equation for H'^^'> with zero energy is "ip^^ = e~^. 

In all cases when one has the factorization relations of the form (1), the super- 
charge operators and the superhamiltonian Hs can be introduced [1]: 
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The operators and H^°^ , H^^^ are called the components of the supercharge and 
of the superhamiltonian, respectively. The relations of the SUSY algebra: 

{Q+,Q-} = Hs (Q+)2 = (g-)2 = [Q^,Hs]=0 (4) 

follow from the definitions (3) and the factorization (1). The last equality in (4) is 
equivalent to the following intertwining relations: 

H^'^S- = S-H^'^ S+H^'^ = (5) 

The eigenfunctions of the Hamiltonians are interconnected by the operators : 

hence their eigenvalues coincide, except for the zero modes of the supercharges. 

Similarly to that, in two dimensions the following pair of Hamiltonians can be 
factorized": 

//(o) ^ -A + {dixf -a3Ax = S+S~ (6) 
//(I) ^ _A + {da? + 2a,d,dsx + cr3idl-dDx = S-S+ (7) 

9, = ^ (i = l,3) A = dl + dl 

This time, the matrix differential operators: 

= - ia2d3 + aid^x + (^zd^X S+ = (S-)^ (8) 

intertwine the Hamiltonians. The superpotential x(x) is an arbitrary real scalar func- 
tion. Note that H^^^ is a diagonal Hamiltonian. The solution of its upper component 
with zero energy is i/jq^^ — e"^. 

We remark that it would be natural to generalize the 1-dimensional factorization 
(1) onto the 3-dimensional models, similarly to (6)-(8). In this case (8) has to con- 
tain a term with 82 multiplied by some Pauli matrix (in order to get a factorizable 
Schroedinger-type Hamiltonian). However, all three Pauli matrices have already been 
used in (8), so a simple direct three-dimensional generalization of (6)-(8) is impossi- 
ble. Some generalizations for 3-dimensional case were discussed in [2], where 4x4 
Dirac matrices were used, and in [3]. But the construction of the last paper is not 
adapted to SUSY-diagonalization of Pauli operator. The 3-dimensional factorization 

"It will be convenient to denote the two-dimensional vector with components R3 as R; Uj, i = 
1, 3 are the Pauli matrices. 
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of the components of superhamiltonian was considered in [4] leading to the Hamilto- 
nian with terms hnear in spatial derivatives (in particular, as a spin-orbital coupling). 
In this case supercharges with 0(3)— rotational invariance were chosen. Here we gen- 
erally do not suppose any spatial symmetry of S^. Another possible 3-dimcnsional 
generalization of (6)-(8) will be discussed in the Section 7 of the present paper. 

Returning back to the factorization relations (6), (7) one could again construct the 
supercharge operators and the superhamiltonian^: 




where 

£(0) ^ _A + (d.xf - Ax H^'^ = -A + (da)' + Ax i = 1, 3 

and is still defined by the first equality of (7). Below we will call H^°^ the diagonal 

component of the superhamiltonian and H^^^ the matrix one. 
From the factorization relations (6), (7) it follows again that 

{Q+,Q-}^Hs; (g+)2 = (g-)2 = [Q^,Hs]^0. (10) 

The last equality of (10) will again be equivalent to the intertwining relations: 

H^^)S- = S-H^'^ S+H^'^ = (11) 

The components of the supercharge operators (8) can be rewritten as: 

where {I, m = 1, 3) 

= Tdi + dix Pi^ = eimQZ en = £33 = €13 = -631 = 1. 

If we plug (12) into the intertwining relations (11) we, after some algebra, will recover 
the usual relations of the two-dimensional SUSY QM [5], [6]: 

Hj^P^n - PrH^'^ P^h\2^H^''^P:^ /,m = l,3. (13) 

^In [5] a different definition was used: tlie positions of tfie components H}^^ and H^^^ in Hg were 
interchanged. For the readers famihar with the two-dimensional SUSY QM [5], [6] we may add that 
this corresponds to choosing a different representation of fermionic creation/annihilation operators. 
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The eigenfunctions of the Hamiltonians H^^\ H^^^ are interconnected by the op- 
erators 5'='=: 

V'«(x; = E-'/'S-^^'\^; E^). (14) 

It follows that the eigenvalues of the Hamiltonians H^^\ if^^^ coincide, except for the 
zero modes of the supercharges. If H'^^^ is not diagonalizablc by constant unitary 
transformations (see Section 3) then the above supersymmetric transformations can 
simplify the derivation of its spectrum by relating it to the spectrum of a diagonal 
Hamiltonian H^^'> [5], [6]. In the most interesting example [7] one identifies H^^^ with 
the Pauli operator [8] of a non-relativistic spin 1/2 particle in external electromagnetic 
field. 

The Pauli operator was studied in the framework of SUSY QM in many papers (see 
[3], [2], [9] - [13]). In these papers the Pauli operator was usually identified with the 
total superhamiltonian, while in the present work the Pauli operator coincides with 
a component of the superhamiltonian (in this sense we continue the approach of 
[7]). An advantage of our approach is the opportunity to consider arbitrary values of 
the gyromagnetic ratio g = — {m contrast to the usual restriction that g = 2). Also, 
we will consider the magnetic fields having non-zero components in the {xi, x^) plane, 
for which the Pauli problem is not diagonalizablc by rotations. We have to mention 
that interesting concrete models of Pauli operator were investigated by means of SUSY 
QM: a particle with spin in the field of monopole [2], and a spin 1/2 particle in the 
magnetic field of a thin straight wire with current [10]. The SUSY-diagonalizability 
of Pauli operators was considered in [12], [13] but only for external fields with definite 
spatial symmetries. 

The paper is organized as follows. In Section 2 we consider the simplest case, 
when the Pauli operator is literally identified with the matrix component of the 
superhamiltonian . The Pauli operator is originally three-dimensional but admits the 
separation of the X2 variable. We derive the expressions for the magnetic field and 
the scalar potential in the Pauli operator and for the superpotential that follow from 
such identification. 

In Section 3 we enlarge the class of the Pauli operators, that can be identified 
with the matrix component of the superhamiltonian, by allowing constant unitary 
rotations, and describe the class of supcrpotcntials for which it is possible. The 
content of Sections 2 and 3 should be considered as a generalization of [7] for the case 
of nonzero electric current. 

In Section 4 we investigate the case of a charged particle in detail. It turns out 
that the separation of X2 variable is not trivial. The class of magnetic field for which 
the Pauli operator of a charged particle can be identified (up to a unitary rotation) 
with H^^^ is described. It is smaller than the one for a neutral particle. 
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In Section 5 we describe a new generalization of the SUSY QM transformations. 
It involves the transition from the usual derivatives to the covariant ones both in 
the components of the supercharge operators and those of the superhamiltonian. We 
also prove the general statement that the quasifactorization relations [5], [6] always 
describe a pair of factorized Hamiltonians, one of which is diagonal. 

In Section 6 we use these generalized transformations to diagonalize the Pauli op- 
erators for more general configurations of magnetic field in the case the gyromagnetic 
ratio is equal to 2. In particular, the component B2 of the magnetic field along the 
X2 direction can be arbitrary. Examples of the Pauli operators allowing the above 
treatment are given in the end of Sections 2, 4 and 6. 

In Section 7 we construct a three-dimensional generalization of the two-dimensio- 
nal factorization relations described in two previous Sections. 

2. Pauli Hamiltonian as a matrix component of the super- 
hamiltonian. 

A non-relativistic particle with spin 1/2 in external electromagnetic field is described 
[8] by the Pauli Hamiltonian'^: 

Hp = -& - lid ■ B{x) + U{x) (15) 
(h^c^2m=l) D = d-ieA{x) 

— * — * — * 

where A{x) is the electromagnetic vector potential; B{x) = i:otA{x) is the external 
magnetic field; U (x) is the external scalar potential (not necessarily an electrostatic 
one); e,fi are the charge and the magnetic moment*^ of a particle. The operator (15) 
acts on the two-component wave functions. 

Let us assume that all external fields in (15) depend on two coordinates (for 
example {xijX^) = x) only. Then the wave function of a particle that moves freely 
along the X2 direction can be written as: 

^(x) = e-'*^"^^(x) (16) 

and the Pauli Hamiltonian, acting on ip{x.), assumes the form: 

Hp = -D^ + {k + eA2(x))2 - ■ fi(x) + f/(x). (17) 

Insofar as the operator (17) is a 2 x 2 matrix differential operator in the two- 
dimensional space {xi, X3), one can try to identify it with the matrix component H^^^ 

"^From this moment on, R denotes a three-dimensional vector with the components Ri, R2, Rs- 
"^We stress that the value of the gyromagnetic ratio g = 2/x/e can be arbitrary. 
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of the superhamiltonian: 

Hp = H^^^ + Eo. 

This idea was set forth in [7], but in that paper it was additionally assumed that 
the magnetic field has no sources: j = 0. In the present paper, on the contrary, we 
consider the most general case of the equivalence of Hp and H'^^\ when the current 
density can be nonzero. Thus one can try to consider a wider class of Pauli operators. 

If such identification is possible, the corresponding Pauli operator can be diagonal- 
ized by the supersymmetric transformations from Section 1, and its spectrum is the 
sum of the spectra of and H^^^ (see (13), (14)). It is possible only if the following 
relations between the external fields in the Pauli operator and the superpotential x(x) 
are satisfied : 

-fxB = {2d,dsx,0,{dl-^,)x) (18) 
C/(x) = (dixT-^k + eA^f + Eo i = l,3. (19) 

Prom the absence of first derivatives in H^^\ and, consequently, in Hp, it follows 
that either e = (a neutral particle) or Ai — As — 0. 

In the equation (19) all terms'^, except {k + eA2)^ and possibly £"0, have no de- 
pendence on k. The k-dependent parts of these two terms cancel only if e = and 
£"0 = k"^, therefore we assume this choice of e and Eq in the remaining part of this 
Section. Prom the Maxwell equation diBi = and (18) one can infer the following 
constraint onto the superpotential x(x): 

ds{3dl - dl)x = 0. (20) 

The general solution of this equation has the form (see Appendix) : 

X(x) = F{x,) + G{-Xi/2 + {^/3/2)xs) + H{-xi/2 - {V3/2)xs) (21) 

where F, G, H are arbitrary thrice-differentiable functions of their arguments. 

The equalities (18), (19) allow us to express all the physical quantities in terms of 
these functions. In particular, 

5i(x) = ^{G" - H") i?3(x) = -i-F" + 1/2(6"' + H")) (22) 

the paper [7], A2 was chosen such that k + eA2 is fc-independent. We will not use this 
idea because A2 is a quantity that enters in the original three-dimensional Pauli operator (15) and 
therefore cannot depend on k. In the following Sections we will discuss cases when x depends on k, 
but for now we assume that it is fc-independent. 
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and the electric current j(x) = rotS(x) has, obviously, only one nonzero component: 

The scalar potential (19) can then also be expressed in terms of F,G,H: 

C/(x) = F'^ + G'^ + H'^ - F'G' - F'H' - G'H'. (23) 

Thus the Pauli operator with the external fields as in (22), (23) coincides with the 
Hamiltonian if*^^) and can be diagonalized by the supersymmetric transformations. 
One can check that in the case j2 — our solution for x(x) coincides with the one 
offered in [7] for the case of a neutral particle. 

Example . Let us choose specific functions F, G, H in (21): 

F{x) = G{x) = H{x) = ax^ 

where a is an arbitrary constant. Then, from (22), 

9a 9a 

-Dl = X3 B3 = Xi. 

The electric current that generates this system of external fields is homogeneous 
and is directed along the axis X2: 

9a 

32 = 



27r/x 

The scalar potential (23) in Hp depends on = x\ + x\ only: 

t/(x) = {d^xf = 1P' 7 = ^ 
and the diagonal component of the superhamiltonian is proportional to a unit matrix: 

since Ax = in this case. Thus the spectrum of the Pauli operator in question 
consists of twofold degenerated levels of the anharmonic oscillator 7p^. 
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3. Unitary rotation. 

In order to enlarge the class of external fields, for which the Pauli Hamiltonian can be 
diagonalized by the supersymmetric factorization, one can identify Hp with the oper- 
ator, which is unitarily equivalent to the matrix component of the superhamiltonian 
[7]: 

Hp = H^^^ + Eo = UH^^^U^ 

where W is a constant 2x2 matrix U — -\- iaa; + — 1; ao, Oij are real 
constants; <? = (cii, (T2, CT3). 

Such rotation obviously doesn't change the spectrum of the matrix operator H^^\ 
but leads to additional freedom because of the presence of the new parameters ao, aj. 
Now the relation between the components of the magnetic field and the superpotential 
X is more complicated compared to (18): 

jiBi = [{al - a'^)Sik + 2{ajao€ijk + akai)]Xk (i, k = 1,2, 3) (24) 

where X = (— 29i93X, 0, (5| — dj)x)- The scalar potential is still defined by (19). 
Taking the divergence of the equation (24), we obtain: 

= iJ^diBi = di[{al - a^)5ik + 2{ajaoeijk + akai)]Xk = 
^-{al + aj-al- al)2dfd3X + (^o - a? - + o;^)(9f - dfds)x + 
+{2aia2 - 2a2ao){didl - dl)x - (2q;2Q;o + 2aia^)2didlx- 

Hence, the superpotential x satisfies the following partial differential equation 
with constant coefficients: 

2{a20iQ — aia^)d\ — (Sckq -\- a\ — 3al — afjd^ds — 2{aia3 + 3a2ao)didl + (25) 



I / 2 2 2 I 2\ o3 



x = o. 



The equation (20) is a partial case of (25) with ao — l;ai — 0. 

It is shown in the Appendix that the supcrpotcntials x, satisfying this equation 
and therefore allowing the identification Hp — H^^^ + Eo, must have the form-'^: 

X = F{tiX3 + Xi) + G{t2X3 + Xi) + H{t3X3 + Xi) (26) 



/If some of ti are complex one can take only such functions F, G, H that the superpotential is a 
real function. 
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where ti,t2, are the roots of the characteristic polynomial, which is determined by 
the operator in (25): 

{2{a2ao — aia3)dl — (Sag + ctf — 3al — al)did3 — 2{aia3 + 3a2ao)didl + 

+2{al -al-al + al)dl) = co{ds - t,d,){ds - t2d,){d3 - tsdi) (27) 

Co = 2(q;q — al — + ccg). 

In the case of a neutral particle (e = but 7^ 0) any choice of ctj and the corre- 
sponding superpotential (26), leads to a Pauli operator determined by (17), (19), (24), 
that is diagonalizable by supersymmetric transformations. 

For a charged particle (e ^ 0) there are stronger restrictions: as we mentioned in 
Section 2, this approach is valid only for Ai — A3 — 0, and, consequently, B2 — 
and ji = j3 = 0. Then, from (24), 

{2a3ao + 2aia2)Bi{x) + {2a3a2 - 2aiao)B3{x) = 0. (28) 

If either of the constants in (28) is nonzero, the Pauh operator can be trivially diag- 

— * 

onalized by a constant unitary rotation, after which B is directed along x^. Hence, 
the interesting situation is when 2a3ao + 2Q;ia2 = 2a3a2 — 2aiaQ — 0. One can check 
that it is possible only in the following two cases: 

1) 0:3 = «! = 

2) a2 = ao ~ 

Unfortunately, in the variant 1) we couldn't find solutions that satisfy (19) with 
the scalar and vector potentials U and A independent on k. 
In the variant 2) the roots ti, ^2, ^3 (27) have the form 

2q;iQ;3 

h = ^ t2^l ^3 = 

0:3 — af 

Hence, taking the reality of the superpotential into account, (26) leads to 

X = F[xi-\ — 5 kX3 + Re giz) z = xi + 1x3 (29) 

\ ai-al J 

where g{z) — 2G{z) is an arbitrary holomorphic function of z. 



4. The dependence of the spectrum on the momentum along 
the axis x^. 

It can be shown that in the case of a charged particle (e 7^ 0) the superpotential (29) 
can satisfy (19) with the scalar and vector potentials U and A = (0, A2, 0) independent 



10 



on k only if F = 0: 

X = Keg{z). (30) 
Then (24) has the following solution^: 

Si = --Im/ S3 = --Re/ A2 = --Re^' + a (31) 
ji II ji 

Let us choose 

2e 

9{z) = go{z) -kz (32) 

4e 

where both qq and Cq do not depend on k. Then from (31) 

£i = --Im^;' 53 = --Re^;' A2 = --Re^; + Co. (33) 

/X 

Note that doesn't depend on k. Now let us prove that U can also be made 
A;- independent. Plugging (30), (32), (33) into (19), we can check that 

U = \g',\' - (^Re^; - Coe)' + E,-(l- - 2Coek. (34) 

By choosing 

Eo^{l- %)k'' + 2CoeA; (35) 
A* 

we can make U /c-independent. Indeed, from (34) 

U=\g',\'-{^j-Reg',-Coe)\ (36) 

Plugging (30), (32) into (6), wc sec that the Pauli operator with external fields 
(33), (36) is equivalent to the diagonal Hamiltonian 



= -A + 



90 - —k 
1^ 



-a^ARego. (37) 



^'A similar expression for the magnetic field was found in [7]. However in the present text we 
allow j2 = diB^ — d^Bi to be nonzero (other two components of the current density are zero because 
of our assumptions). This is possible if go{z) has singularities, i.e., it is not an entire function (see 
Example below). Another difference of the approach of [7] is that it did not allow to obtain the 
dependence of energy levels on k, in contrast to ours. 
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Note that k enters into H^^^ nontrivially. The last term in H^^^ is zero everywhere 
except the points where g has singularities. 

Example . Let us choose a specific function go{z) in (32): 

go{z) = alnz 

a being a real constant. Then we readily obtain from the corresponding general 
expressions that 

A _ ^axi D _ 4" ^1^3 n _2axl-xl 
+ Uo -fc>i — J— -D3 



j2 = did{xi)S{x3). 

A* 

This electric current density can be interpreted as a "dipole current" - two thin 
linear currents in converse directions, the distance between them being small. 
In this case, (36) turns into: 



P^ \ pp-^ ) 



Prom (37) it follows that the spectrum of our Pauli operator coincides with the 
spectrum of the diagonal Hamiltonian H^^^ ^ which has the form: 



a 2ek 



2 



2T:a5{xi)5{x^). 



5. More general factorization method. 

In the Pauli operator (17) each derivative enters in the covariant derivative combina- 
tion Dj = dj — ieAj. Therefore, it seems natural to use the covariant derivatives in 
the operators (6)-(8) too'*: 

S^ = TDi-ia2Ds + a,d3X + (T3diX S+ = iS-)l (38) 

The so defined intertwining operators generate the following pair of factorized 
Hamiltonians: 

i/(o) = S+S- = -T>' + {day-ea2B2-asAx (39) 
i^W = s-S+^-T>^+{day + ea2B2 + 2a,did3X + 'J3{d!-di)x (40) 
B2 = d^A.-diAs i = l,3. 



''Note that in the simplest particular case x = the operators (38) can be reduced to the SUSY 
generators given in [3], [9] by the multiplication onto ai or 0-3 and a redefinition of axes. 
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The Hamiltonians (7) obey the intertwining equations similar to (11): 

and their eigenfunctions are interconnected. 

Similarly to Sections 2,3 one can try to identify H^^^ with a Pauli operator up to 
constant unitary rotations: 

Hp = H^^^ + Eo H^^^ = UH^^^U+ (41) 

and H^^^ is unitary equivalent to a diagonal Hamiltonian: 

= V+^(°)V (42) 

If such identification is possible, the corresponding Pauli operator can be diagonalized 
by the above supersymmetric-like transformations and its spectrum is the sum of the 

spectra of "* and "* Below we shall show that the more general factorization 
method presented in this Section allows us to diagonalize a broader class of Pauli 
operators than the usual supersymmetric methods described in the previous Sections. 
From (39)- (42) one may infer the new factorization relations: 

^(0) = S+S- H^^^ = S'S^ S- = US'V^ S+ = VS+U-^ = (^5")^ (44) 

Note that from (44) we can construct the SUSY algebra relations, similar to (3), (4), 
or (9), (10), and thus obtain a new realization of SUSY QM in two dimensions, that 
differs from the one given in [5], [6]. 

Prom (44) the following intertwining relations can be derived: 

and the eigenfunctions of H^^\ H^^^ are interconnected by the operators 5*^: 

V^W(x;K)=£^-'/'^-V'^°Hx;K). 

Thus the eigenvalues of the Hamiltonians H^^\ H^^^ again coincide up to the zero 
modes of the supercharge components S^. 

Now let us parameterize the matrices S"^ in the following way, which is similar to 
(12): 

I I) ^^=(1 I)- («) 
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Plugging these operators into (44), (43), one can infer the following quasifactorization 
relations: 

S^'^=QtQi ^'^^P^Pr- h12 = QTQI + Pi-p:^. (46) 

Pi^Qi =QiPi- = ^ /,m=l,3. 
Prom them we can deduce an analog of the intertwining relations (13): 

H^^P^ = ^rl^°^ Pi'Hi:^-i'^P^ i,m = l,3. (47) 
One may also note that H^^^ can be decomposed as: 

hS1-H+^1 M^l^QTQl. M^'l^Pi-P;^ (48) 

such that H^JiM!^,.], — Mzim—mn — 0) ^^'^ relations (47) can be rewritten in the 
form: 

M12p^ - PfH^'^ Pi'3!'1 = H^'^P:: Z,m = l,3. (49) 

The relations (46)- (49) have exactly the same form as the well known SUSY QM 
quasifactorization relations and the SUSY QM intertwining relations [5], [6]. The 
former are reduced to the latter in case U — V — I and Ai—Ay, — ^. 

From the above algebra we see that the any quasifactorization relations in fact 
describe a pair of factorized (in the sense of (44)) Hamiltonians, one of which is 
diagonal. This is a rather general statement, because the relations (46)-(49) follow 
from the definitions (43)-(45), no matter what form the components of the supercharge 
Q±,P± have. 

6. The cases when if^^^ is diagonahzable. 

Eq. (41) is equivalent to the following pair of equations: 

= [{al - d'^)5ik + 2{ajaoeijk + akai)]Yk (i,/c = l,2,3) (50) 
C/(x) = (dx)' - (k + eA^)' + Eo (51) 

where Y = {-2did^x, -eB2, -{dj - dl)x)- 
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Let us consider for simplicity the case when the gyromagnetic ratio 2///e is equal 
to^ g = 2 and ao = 02 = 0. Then (50) is reduced to (24), which is solved in the 
Section 3, the solution being given by (26). Assume further that =1; ai — 0, i.e., 
U = a^. Then, Eq. (26) is reduced to a partial case of (29): 

X ^ F{xi) +Reg{z) z ^ Xi + ixs. (52) 

Plugging the superpotential (52) into (50) we obtain: 

2 2 1 2 1 

Bi = — Im g" B3 = — Re g" - -F" A2 = — Re g' - -F' + const, 
e e e e e 

The first component H^^^ of the superhamiltonian is unitarily equivalent to a diagonal 
Hamiltonian in the following two-' cases: 

a) there is some constant Cb such that B2 = C^Ax 

b) Ax = 0. 

Unfortunately, in the case a) it is impossible to satisfy (51) when both U and A2 
are /c-independent. In the rest of this Section we will consider the case b). Then 
the Hamiltonian is already diagonal; therefore, in the rest of the article we will 
assume that V = I. 

If the superpotential x (52) satisfies the condition Ax = 0, it follows that F{xi) 
is a linear function. One can check that we can set F{xi) — without narrowing the 
class of Pauli operators that can be diagonalized by our factorization method. So, 

X^Regiz). (53) 

This is the same superpotential as in (30). However, this time g{z) cannot contain 
singularities anywhere except at the infinity, because otherwise Ax would be nonzero. 

The magnetic field components -81,-83 are obtained from the same equation as 
(24), therefore they will have the form*^ (31): 

Si = --Im/ B3 = --Reg" = --Re^' + C. 

e e e 

The component B2 can be arbitrary. Similarly to the Section 4, we choose the super- 
potential and the constant in the magnetic potential in the form (32): 

g{z) = go{z) - 2kz (54) 



'This is the first instance when we fix the gyromagnetic ratio to the same value g = 2 that 

appeared in [3], [9]. However, in this text we consider a much larger class of the configurations 
of external fields compared to [3], [9], because, in contrast to them, we allow the magnetic field 
component in the (a: 1,2:3) plane to be nonzero. 

•' We do not consider the case when Ai = As = that has already been investigated above. 

'^We keep taking into account that the gyromagnetic ratio 2/x/e is fixed at 2. 
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Then the magnetic field components Bi, B3 take the form (33): 

B,^--Img'^ Bs^--Reg'^ (55) 
e e 

while B2 remains arbitrary. We get the scalar field U{x.) from (51), which coincides 
with (19). If we assume that Eq has the form (35): 

Eq = -3k^ + 2Coek 

then the scalar field will be defined by (36): 

U=\gX-(2Reg'o-Coey 

Thus, same as in the Section 4, all physical fields in the Pauh operator are k- 
independent. 

Plugging (53), (54) into (39), we see that the spectrum of the Pauli operator coin- 
cides with the spectrum of the Hamiltonian 

//(o) ^ _d2 _ ^^^^ ^ _ 2^|2_ 

After a unitary rotation by V = (a2 + cr3)/\/2, H^^^ will turn into a diagonal Hamil- 
tonian analogous to (37): 

^(0) ^ _ ^^^^ ^ _ 2fc|2_ (5g) 

Example . Let us choose a specific function g{z) in (53): 

/ ^ 3 

9i^) = • 

Then we readily obtain from the corresponding general expressions (55)-(56) that 

- — (X3,0,xi) + (0,^2,0) j = 6 
e 

U = - {u{xl - xl) + canstf H^^^ = -D^ - (73^2 + 

7. A three-dimensional factorization. 

In this Section we will show that the above factorization method allows a three- 
dimensional generalization. It is realized by the operators 

R^ = Tff-D + Vo{x); R+ = {R-y (57) 



2k 
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that can be used for the factorization' of the following pair of Hamiltonians: 

= R+R~ ^ -D'^ -a ■ (dVo + eB) + V^; (58) 
= R-R+ = -D^ + a-{dVo-eB) + V^, (59) 

where again D = d — ieA; B = rot A; Vo{x) is an arbitrary real scalar function. 
The operators (58), (59) can then be treated in the same way as in the Section 5: 
can be made unitarily equivalent to a Pauli operator in three dimensions (15): 

~(1) ^ ^ ru-^ 

Hp^H +Eo H ^UH^^^U+, (60) 
and ' to a diagonal Hamiltonian H : 

^(0) = V+H°^V. (61) 

Now let us prove that the operators (57) - (59) are indeed a generalization of (38) 
- (40). Namely, assume that Vq and A depend on two coordinates (xi, X3) only. Then, 
similarly to Section 2, the wave function of a particle along the X2 direction has the 
form (16), and the operators (57) - (59), acting on ip{x), assume the form: 



+ Vo; (62) 



aiDi + asDs — i(k + eA-2)a2 

^(1)^^(1) = -B^ + {k + eA2f + a-{dVo~eB) + Vo^ = R-R+. (64) 

The Pauli operator and the diagonal Hamiltonian in (60), (61) will also transform to 
two-dimensional operators. 

Now, if we assume that Vq — d^x and A2 — ^(<9iX — k) then one can check that 
(62) - (64) can be rewritten in terms of the operators from the Section 5: 

^+ = ai5+; R- = S-a^; ^(0) = ifW; H^'^ = a^H^^^a^ 

where are defined in (38), and H^^\^ H^^^ in (39), (40). Then one can check that 
the two-dimensional versions of (60), (61) arc equivalent to (41), (42), if we set'": 

U = Uai] V = V. 

Thus, we have proved that the operators (57) - (59) are indeed a generalization 
of (38) - (40). Such generalization may allow one to diagonalize a more general class 
of three-dimensional Pauli operators. 

'The factorization (57) - (59) was proposed in [14] for the study of massless Dirac operators in 
Euclidean spacetime. 

'"Note that the choice of the unitary matrices V = /; U = <T2, that we have widely used in Sections 
3-6, is thus equivalent to the simplest choice in the new formalism of this Section: V = U = I. 
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Appendix. 

It is straightforward that any differential operator of the form 

codi + cidld^ + C2dsd^^ + c^dl 

with Co 7^ can be factorized into a product of differential operators of first order: 

co{dz - tidi){d3 - t2di){d3 - Udi) 

where ^1,^2,^3 are the roots of the polynomial CqX^ + Cix"^ + C2X + C3. Assume for 
simplicity that these roots are different. Then the following statement is true: if 

(^3 - hd^){d3 - t2di){d3 - hd^)f = (65) 

then / = F{tiX3 + xi) + G{t2X3 + xi) + H{t3X3 + xi). 

Proof : if (65) is satisfied, then there exists a function fo{tiX3 + Xi) such that 

{d3-t2d,){d3-t3di)f = fo{tlX3 + Xi). (66) 

Insofar as ^1,^2,^3 are assumed to be different, one can show: there exists a function 
F{tiX3 + xi) (see below) such that (66) is equivalent to the following equation: 

(^3 - t2di){d3 - t3d,){f - F{tiX3 + xi)) = (67) 

where (^3 - t2<9i)(<93 - t3di)F{tiX3 + xi) = /o(tiX3 + Xi). 

Similarly, from (67) it follows that there exists a function ^(^22^3 + Xi) such that 

(^3 - t3di) (/ - F{hX3 + Xi) - G{t2X3 + Xi)) = 0. 

Therefore, there is H{t3X3 + Xi), such that 

/ - F{tiX3 + xi) - G{t2X3 + xi) - H{t3X3 + a;i) = 0. 

Q.E.D. 

Note that in case when t2 — ^2 for / to be real it is necessary that G = H, i.e. 
f = F + 2ReG. 
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